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ABSTRACT 

Dark-matter halos grown in cosmological simulations appear to have central NFW-like density cusps 
with mean values of dlogp/dlogr « — 1, and some dispersion, which is generally parametrized by the 
varying index a in the Einasto density profile fitting function. Non-universality in profile shapes is 
also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non-universality, 
at any given mass scale, is an intrinsic property of DARKexp, a theoretically derived model for 
collisionless self-gravitating systems. We demonstrate that DARKexp—which has only one shape 
parameter, cf>o- fits the dispersion in profile shapes of massive simulated halos as well as observed 
clusters very well. DARKexp also allows for cored dark-matter profiles, such as those found for dwarf 
spheroidal galaxies. We provide approximate analytical relations between DARKexp </)$, Einasto a, 
or the central logarithmic slope in the Dehnen-Tremaine analytical 7-models. The range in halo 
parameters reflects a substantial variation in the binding energies per unit mass of dark-matter halos. 
Subject headings: dark matter — galaxies: halos — galaxies: kinematics and dynamics 


1. INTRODUCTION 

The basic shape of N-body simulated cold dark matter 
ha los is now w el l est ablished. The NFW density profile 
(INavarro et all 119971) is a remarkably good representa- 


tion (e. 


2005, 


., .Tasitsiomi et al.l 12004 Dicmand ct al. 

INavarro et al.1 l2010fc iDiemand fc Moore 


2004, 


2011 


Ludlow et al.l I2011D . especially given its simplicity: it 


has no shape parameters. However, a closer examination 
of the profile shapes of relaxed halos shows that NFW 
does not fully capture the density run. Most impor¬ 
tantly, while the NFW profile has a central density cusp, 
with logarithmic slope d log p(r)/<i log r = —1, recent 
high resolution simulations have shown that the central 
cusp tends to beco me somewhat sh allower at smaller and 
smaller radii (e.g., INavarro et al.ll20 1f)D. prompting the 
use o f an Einasto oemodel (jNavarro et al.ll2004 lEinastol 
I1965D . The Einasto profile fits halos better at most radii, 
not just in the center, suggesting that the additional pa- 
rameter a is a real characteristic of d ark matter halos 
dMerritt et al.ll2 0()fil: iKlvnin et al.ll2014l ). 

iLudlow ct al.1 (12013 ) find that the mean value of a 
ranges from about 0.17 for M200 = 10 11 /i _1 M 0 ha¬ 
los to 0.2 1 at M?nn = 10 15 h~ x M 0 , which is consis¬ 
tent with I Cao et ahl (2008) who find a relation between 
Einasto’s a and the ratio of the linear density thresh¬ 
old for collapse to the rms linear density fluctuation at a 
given mass scale. In addition to the mass-dependent a, 
there is significant scatter of around 0.05 in a. Moreover, 
halos with a in the range 0.05 to 0.45 are also realize d 
in simulations ((Ludlow et al.1 [2013 : iKlvpin et al.ll2014l ). 
These works show that the shape of the density profile, 
including the central logarithmic slope, is tied to the ini- 
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tial conditio ns and evolutionary history of the halos in 
simulations. ILudlow et al.l (|2013l) suggest that the mass 
profiles of halos reflect their accretion history. 

The spread in a inde x is observed n o t just for main 
halos, but also subhalos. iSnringel ct al.1 ( 2008 ) find that 
the profiles of subhalos do not converge to a central 
power law but exhibit curving shapes well described 
by the Einasto form with a in the range 0. 15 to 0.28. 
This i s consistent with the findings of IVera-Ciro et al.l 
(I2013D , who note that simulated satellite dark matter ha¬ 
los around galaxies are better fit wit h Einasto profiles of 
index a = 0.2 — 0.5. iDi Cintio et al.l (|2013t) find a corre¬ 
lation between a -1 (ranging from 0.1 to 1) and satellite 
dark matter halo mass, i.e., with higher masses corre¬ 
sponding to smaller values of a, contrary to the trend 
found for more massive halos. The diversity in the den¬ 
sity profiles of these satellites is attributed to tidal strip¬ 
ping, which acts mainly in the outer parts of the sub¬ 
haloes. 

The conclusion is that simulated dark matter halos 
do not form a strictly self-similar family of systems as 
was originally suggested, but display non-universality 
dMerritt et alJM . 

Neither NFW nor Einasto have any theoretical un¬ 
derpinning s; they are purely emp i rical. In a series 
of papers dHiorth fc Willi amsl 1201(1 : iWilliams fc Hiorthl 
1201(1 : Williams et al.l 1201(1 l2nTl ~ we have proposed a 
theoretically based model, DARKexp, which has been 
shown to allo w good fits to a few simulated halos 
dWilliams et al.l 120101 I2014D and observed galaxy c lus¬ 
ters dBeraldo e Silva et al.ll2013l : lUmetsu et alJl2dl5l) . 

DARKexp has one shape parameter, and its density 
profiles are very similar to Einasto within a limited range 
of shape parameters. Taking the broader range of achiev¬ 
able values into account, DARKexp suggests that there 
should be a dispersion in the shapes of the density profiles 
between systems, parametrized by the shape parameter. 
At small radii, but sufficiently above the resolution limit 
achieved in cosmological simulations, these density pro¬ 
files differ significantly. Some have nearly flat cores, while 
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others are steep, d log p(r)/d log r « —2. In other words, 
unless halos have identical normalized central potentials, 
DARKcxp predicts that the density profiles should not 
be universal, and that a logarithmic slope of —1 at small 
radii robustly resolved in simulations is not a univer¬ 
sal attractor of halos with different formation histories. 
Testing this would be important because DARKexp is 
not a fitting function but is derived from basic physics. 

We summarize the rationale and properties of DARK¬ 
exp in Section [2j While DARKexp has the advantage 
of being derived from first principles, the resulting den¬ 
sity profiles unfortunately cannot be expressed analyti¬ 
cally. We therefore also relate the properties of DARK¬ 
exp to those of phenomenological fitting functions, i.e., 
generalized Hernquist, Dehnen Tremaine, and Einasto 
models. We discuss non-universality of simulated and 
observed massive dark-matter halos in Section 3 and 
show that, in cases where exact fits are not required, 
DARKexp, Einasto, or Dehnen Tremaine profiles pro¬ 
vide equally good representations, with considerable dis¬ 
persion in their respective shape parameters. We relate 
observed, simulated and theoretical properties of lower 
mass halos in Section 4 and conclude in Section 5 

2 . DARKexp AND THE CENTRAL 
POTENTIAL-DENSITY SLOPE RELATION 

DARKexp is a theoretical model for the differential 
energy distribution in a self-gravitating collisionless sys¬ 
tem. Specifically, defining N(e) = dM/de , where e is 
the normalized absolute particle energy per unit mass, a 
statistical mechanical analysis leads to 

N(e) oc exp(0o - e) - 1 (1) 

(IHiorth & Williams! 1201 0(1. The central dimensionless 
potential, 0 o, is the only free shape parameter of the 
model. Some of the global properties of the model are 
analytic (see Appendix A), while the phase space distri¬ 
bution function and derived properties, such as the den¬ 
sity profile, must be obtained numeric ally using an iter¬ 
ative procedure dWilliams et all [2010 ,). Apart from the 
assumption of isotropy (which leads to spherical struc¬ 
tures), there are no additional assumptions involved in 
this computation. The scale parameters of the corre¬ 
sponding density profile follow from the resulting profiles 
themselves. DARKexp is a maximum entropy theory and 
as such describes only the final, equilibrium states of sys¬ 
tems, and not the evolutionary paths that brought them 
to the final state. The density profiles apply to isolated 
systems, i.e., there is no cosmological context and so no 
concept of a virial radius or formation redshift. 

While the asy mptotic logarithmic dens ity slope is 
— 1 for r —t 0 (IHiorth fc Williams! I2010f b in reality 
the central logarithmic slope robustly resolved in sim¬ 
ulations depends on 0o: 0o < 4 systems have in¬ 
ner logarithmic slopes shallower than — 1 , while 0 o > 
5 have inner logarithmic slopes between —1 and —2 
((Williams fc Hiorthl 12010( 1. DARKexp was not designed 
to fit the properties of dark-matter structures. Never¬ 
theless, the energy distributions resulting from this first 
principles derivation are excellent representations of nu¬ 
merical simulations and the density profiles are remi- 
niscent of Einasto a profiles for values of 0n around 4 
dWilliams et alJfeOlOH . Note that in the IWilliams et all 


(l2014f) formalism, the DARKexp density run is inde¬ 
pendent of any velocity anisotropy. DARKexp also 
provides very good f its to observed galaxy c lusters 
ilBcraldo e Silva e t al.| (201 31: lllmetsu et all 12015 b and 
globular clusters (IWilliams et al.ll20l2 1. 

We here compare DARKexp to phenomenological fit¬ 
ting functions, notably the generalized Hernquist, the 
special case of the Dehnen-Tremaine 7 —model, and the 
Einasto profile (see Appendix B for details) 0 We pro¬ 
vide the correspondence between the shape parameters 
of these analytic fitting functions and the dimensionless 
parameter of DARKexp, 0 O . These results can be used 
in any situation where fitting with DARKexp is required. 

By construction, the DARKexp energy distribution 
is discontinuous at the escape energy, and as a c onse¬ 
quence the density falls off a s r -4 at large radii (|.Taffel 
Il987t IHiorth fc Madsenlll99ll ). When searching for phe¬ 
nomenological or analytical profiles that may be useful 
appro ximations t o DARKexp, we recall that the gener¬ 
alized iHernauisH (1199011 model has the same outer loga¬ 
rithmic slope of —4 as DARKexp, and also allows for a 
range of central logarithmic slopes, 

p oc r -7 (l + r ^)(i' _4 )/^ ) ( 2 ) 

where /3 governs the sharpness of the transition between 
the asymptotic inner logarithmic slope of —7 and the 
outer logarithmic slope of —4. For /3 = 1 the resulting 
7 -models, 

p oc r -7 (l + r) 7-4 , (3) 

are l argely analytical (iDehnenl 119931 : I Tremaine et al.l 
Fl99l . 

Another phenomenological model, which provides an 
excelle nt fit to numerica l simulations, is the Einasto 
model dNavarro et al.ll200ll : it is controlled by the shape 
parameter a, 


p oc exp 


— (r“-l) 

a 


(4) 


As shown b y IWilliams fe Hiorthl (120101 1 and 
IWilliams et al.l ( 2010lf . DARKexp halos have a range 
of logarithmic slopes at small radii robustly resolved in 
simulations depending on 00 ; models in a range around 
00 ~ 4 are well represented by Einasto models with a in 
the range 0.15-0.2. 

To investigate the relations between the shape param¬ 
eters describing DARKexp and phenomenological mod¬ 
els we now turn to a detailed comparison of their density 
profiles (see Appendix B). For DARKexp of a given 0o we 
find an Einasto model or 7 -model that matches DARK¬ 
exp (i.e., has the same value as DARKexp) in either den¬ 
sity or log log density slope. The anchor radius at which 
the matching is done is chosen to be log(r a /r_ 2 ) = —1.5, 
which is the t ypical resolutio n limit in cosmological sim¬ 
ulations ('e.g.- lNavarro et aT1l2010D . We choose to anchor 
the profiles at a fixed radius rather than fitting the (log) 
density or slope as a function of (log) radius, over some 
range, for visual clarity. 

Figure [T] summarizes the results of the matching be¬ 
tween DARKexp and phenomenological models. It shows 


4 For convenience we are making numerically generated models 
for <j >0 = 1-10 available for download at the DARKexp website 
http: //www.dark-cosmology.dk/DARKexp/ 
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Figure 1. relation (top) and </>o - 7 relation (bottom) for 

DARKexp. The solid curves are for match in density, the dashed 
curve are for match in logarithmic slope at log(r a /r— 2 ) = —1.5. 
The dotted curves represent analytical approximation to these rela¬ 
tions (Equations [3] and [H]). In the top panel, the black symbols rep¬ 
resent simulated halos (Section 3.1, Figure[3}. In th e bottom panel, 
the black symbols are clusters of galaxies from INewman et al.l 
12013al l. obtained as discussed in the text (Section 3.2). 

the relations between 00 and either 7 or a, matched in 
either density or logarithmic slope. We also provide an¬ 
alytical approximations to these relations. For the 7 - 
rnodel 

7 ss 3 log - 0.65; 1.7 < 0 O < 6 , (5) 

which suggests 0o ~ 3.5 for 7 = 1. For the Einasto model 
a « 1.8 exp(— 0 o/ 1 . 6 ) ( 6 ) 

and 0o ~ 3.8 for a = 0.17. 

Comparisons between DARKexp and the empirical 
models, using the above analytical relations, are shown in 
Figure [2] We will return to the central potential-density 
slope relation in the following sections. 

3. NON-UNIVERSALITY OF MASSIVE HALOS 
3.1. Simulated dark matter halos 

Until about 5-10 years ago simulation results were con¬ 
sistent with halos being universal in shape. Now, there 
is considerable evidence to the contrary. How does the 
DARKexp expectation of non-universality compare to 
simulations? We here analyze recent simulations, show 
that DARKexp fits the density profiles at least as well as 
Einasto, and quantify the non-universality. 

We use cluster halos identified in the Bolshoi simula- 


Figure 2. Comparison of DARKexp with Einasto (top) and 7 - 
models (bottom), using the analytical <j>0~ a relation (Equation [71 
top panel) or the analytical <j>0 — i relation (Equation [5] bottom 
panel). DARKexp models with <j >0 = 2, 4, 6 (from bottom to top at 
small radii) are shown in blue (solid curves: density profile; dashed 
curves: logarithmic density slope). 

tior@, a dark matter-only simulation run in a 250h -1 Mpc 
box (WMAP5 cosmological parameters) with 2048 3 par¬ 
ticles, which corresponds to a mass reso l ution of 1 .35 x 
10 8 /i _1 Mq (for more details, see lKlvpin et al.ll2011[ ). For 
the purpose of our analysis, we select well-resolved equi¬ 
librium halos with at least 5 x 10 5 particles within the 
virial radius r v (defined as the radius of a sphere, whose 
enclosed density is 360 times larger then the background 
density). The state of equilibrium is assessed using two 
diagnostics: the virial ratio and the relative offset be¬ 
tween the minimum of the potential and the halo mass 
center. We find 261 equilibrium halos defined by the 
virial ratio less than 1.3 and the mass center offset less 
than 0.05t" v , which are two com monly used relaxation 
criteria proposed bv iNeto et al.1 (120071) . This selection 
is done in order to make a meaningful comparison with 
DARKexp, which by definition only applies to equilib¬ 
rium systems. 

We calculate the density profiles in 30 spherical shells 
equally spaced in a logarithmic scale of radius spanning 
the range from r m i n to r v , w here r m j n = 0.015r v is the 
convergence radius defined bv lPower et al.l (2003 ) and es¬ 
timated for the minimum number of particles per halo. 
The halo centers are given by the minimum of the poten- 

5 The simulation is publicly available t hrough the MuItiD ark 
database (http://www.multidark.org). See lRiebe et al.l 1 2013( 1 for 
details of the database. 
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Figure 3. Spherically averaged density profiles of cluster-sized equilibrium dark matter halos formed in cosmological simulations. The 
profiles are split into three subsamples (columns) corresponding to three ranges of the best-fit shape parameter a of the 7 profile (top), 
Einasto profile (middle) or the best-fit central potential </>0 of DARKexp (bottom). The black curves are the 7 (top), Einasto (middle), or 
DARKexp (bottom) profiles for the mean 7 , a or <j)o in the corresponding halo sample (solid curves: 7 > 1.44, a < 0.164, (j )0 > 4.1; dashed 
curves: 1.30 < 7 < 1.44, 0.164 < a < 0.210, 3.6 < (j )0 < 4.1; dotted curves: 7 < 1.30, a < 0.210, $0 < 3.6). The lower panels in either row 
show the residuals from the best fitting profiles. DARKexp provides as good fits to the simulated halos as the Einasto profile and accounts 
for the range in properties revealed here. 


tial. DARKexp, Einasto, and Dehnen-Tremaine models 
are fitted to the density profiles of the halos by minimiz¬ 
ing the following function 

^shell 

X = [In Pi - Inpmodei(ri)] 2 , (7) 

i=l 

where n s heii is the number of shells and ?y is the mean 
radius of pa r ticles inside the *-th shell. Following 
iLudlow et al.l (|2011h . we restrict all fits to the radial 
range r min < r < 3r_2, which prevents including dy¬ 
namically less relaxed parts of the halos. Here, r _2 is 
defined as the radius at which dlog p{r)/d log r = —2, 
and P —2 = p(r_ 2 ). We first estimate r _2 by fitting an 
Einasto profile in all 30 shells, i.e., within the radial range 


r m i n < r < r v . This sets the restricted fitting range to 
be the same in both cases. Next, r_2, P-2, and 7, a or 
(f >0 are obtained from a restricted fit to either model. In 
fitting DARKexp, we used density profiles computed on 
a grid of </>o, ranging from 2.0 to 8.0 with an interval of 
0 . 01 . 

Figure [3] shows the density profiles of all halos selected 
for the analysis, in three ranges of the shape parameter, 
i.e., a for the Einasto model, 7 for the Dehnen-Tremaine 
model, and <po for DARKexp. The ranges were chosen 
such that there are equal numbers of halos in each of 
the three bins. The mean value of 7 is 1.36 with an rms 
scatter of 0.18 while the mean value of a is 0.19 (rms 
scatter 0.05). For the mean value is 3.9 (rms scat¬ 
ter 0.6). There is no significant correlation in these pa- 
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rameters with the mass over the narrow range of masses 
considered (10 14 - 10 14 8 M 0 ). The lack of any difference 
between residuals from the best-fitting Einasto, Dehnen- 
Tremaine, and DARKexp profile leads to the conclusion 
that these models recover spherically averaged density 
profiles of simulated cosmological halos equally well. 

The second important conclusion is that, in the narrow 
mass range considered, there is a range of density profile 
shapes among simulated systems, and that the DARK¬ 
exp family of models captures that spread very well. Fig¬ 
ure Q] (top) illustrates this pictorially; it shows the best- 
fit shape parameters of the Einasto and DARKexp fits 
to the simulated halos; fo and ol range between 2.5-6 
and 0.05-0.35, respectively (7 ranges from 0.8 to 1.7). 
The scatter in a is about 0.05. As expected (Section 2), 
the data show a tight relation between the central poten¬ 
tial f 0 and the shape parameter a of the Einasto profile, 
which traces the a-fo relations obtained from matching 
Einasto and DARKexp density profiles. 

3.2. Observed clusters of galaxies 

iLimousin et al.l (2008) have emphasized the wide range 
of central logarithmic slopes inferred observationally for 
different clusters of galaxies. DARKexp applies to equi¬ 
librium collisionless systems, so in addition to dark mat¬ 
ter halos it can also be applied to galaxies and clusters 
that host baryons, regardless of the physical processes 
that led to the formation of the final equilibrium sys¬ 
tems, provided that the evolution and the present state 
of these systems are not dominated by two-body en¬ 
counters, and the systems are not far from being spher¬ 
ical. We here show that relaxed galaxy clusters exhibit 
non-universality that can be described by the DARKexp 
model. _ 

iNewman et al.l (2013a) constructed spherically aver¬ 
aged density profiles of seven massive relaxed clusters, 
using central galaxy kinematics on the smallest scales, 
strong lensing on intermediate scales and weak lensing 
on large scales. To separate out dark matter density pro¬ 
files they subtracted the stellar contribution of the cen¬ 
tral galaxy assuming constant mass-to-light ratio for the 
stellar population. The remaining dark matter profiles 
were fit with cored NFW and generalized NFW (gNFW) 
forms. 

The 7 -model we use here and the gNFW differ in the 
asymptotic outer logarithmic slope (—4 vs. —3), so the 
differences are confined mostly to radii outside r_ 2 . At 
r < r _2 the two models are parametrized similarly, there¬ 
fore one can compare th e two models, if one lea ves out 
the outermost radii. The INewman et al.l (l2013af) cluster 
profiles extend about 2 decades in radius interior to r_ 2 , 
and about 0.75 of a decade exterior to r_ 2 , so the bulk 
of the radial r ange is usable in the co mparison. 

We took the INewman et al. C2013a.) inner logarithmic 
slope for each cluster and set it equal to 7 . This sup¬ 
plies the vertical axis coordinate in Figure [D (b ottom ). 
We then fit DARKexp directly to INewman et al] (l2013af) 
gNFW fits, excluding the outer radii beyond which the 
logarithmic slope is steeper than —2.7 (the fits are in¬ 
sensitive to the truncation logarithmic slope for values 
between —2 and —2.8). This gave us fo, plotted along 
the horizontal axis. The seven galaxy clusters are plotted 
as solid points in Figure [T] (bottom). They fit well the 
fo ~7 relation matched in slope and are consistent with 


our proposed analytical approximation (the dotted red 
line). 

From Figure [T| we conclude that clusters are fit with 
a range of <^o’s, suggesting non-universality. Observed 
relaxed clusters of galaxies appear to follow DARKexp; 
and seem to prefer fo values about 2-3.5, on average 
smaller than cluster-sized dark matter halos in pure dark 
matter simulations, which are fitted with values of fo in 
the range 2 - 6 . 


4. DWARF GALAXIES 

So far we discussed dispersion in density profiles among 
massive pure dark matter halos. We next consider the 
differences between pure dark matter halos and observed 
low mass systems. 

The central density profiles of observed galaxy clusters 
appear to be, on average, shallower than those of clusters 
in pure dark matter simulations. Likewise, the central 
logarithmic slopes of observed dwarf galaxies are shal¬ 
lower than simulated dark-matter subhalos. In this sec¬ 
tion we discuss the relative energies between the steeper 
pure dark matter and shallower observed structures. We 
will show that the differences amount to a large fraction 
of the total binding energies of the haloes, demonstrat¬ 
ing that the differences in central logarithmic slopes or 
fo represent significant differences between the halos in 
terms of global energetics. Moreover, assuming the dif¬ 
ference between two states or halos is brought about by 
some dynamical or baryonic proces, we compute what it 
would take to induce a transformation between the two, 
such as the cusp-core transition. Calculating relative 
(fractional) energetics allows us to bypass the uncertain¬ 
ties and unknowns associated with the specific physical 
processes that would bring about the transformation. 


4.1. Lack of cusps in dwarf galaxies 


The differences between the NFW and Einasto de¬ 
scriptions of pure dark matter halos are small in 
comparison to the differences between dark mat¬ 
ter profiles of subhalos/satellite halos and observed 
dwarf galaxies which may have central cores. Most 
dwarf galaxies appear to exhib i t sign i ficant c entral 
density cores dFlores &; Frirriackl 1 1 99-Ii iMoorel 11994 
Strigari et al.l 120061: iGilmore et al.|_ 2007t Walker et al 


20091 : iBoylan-Kolchin et a l . | 20111: ^Mker^^Penarrubia 

lAmorisco fc Evans! 120121: lA morisco et al.l 120131 : 


2011 


.Tardel fc Gebhardtl 120131 [Collins etal.ll2014l f. Whet her 
all dwarfs have a universal profile or not is still unclear 
because of the large uncertainties associated with the 
meas urement of the cent ral logarithmic density slopes 
(e.g., iStrigari et al.l [2014 11: conclusions seem to depend 
on the samp l e and the methods used. For example, 
lAdams et al.l (: 2014 11 claim that dwarfs are consistent 
with having a universal central logarith mic density sl ope 
of —0. 63, with a scatter of 0.28, while lJardel fc Gebhardtl 
(120131) claim that individual dwarf spheroidals show con¬ 
siderable scatter around an average logarithmic slope of 
- 1 . 


For reviews of th is s mall scale controvers y, _ see 

IWeinberg et al.1 (|2013lf and lPontzen fc Governatol (1201411 . 
The suggested solutions to this set of problems are ei¬ 
ther baryonic or particle in nature. War m dark mat- 
ter does not seem to be the solution: iMaccio et al.l 
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(l20Hl20l3h and iShao ct all (|2013T ) find that warm dark 
matter leads to cusps just like cold dark matter does, 
although with very small central cores which however 
are astrophysically uninteresting and much smaller than 
what is cl aimed to be p resent in dwarf satellites. In¬ 
deed, iLovcll ct~. i (12014ft conclude that warm dark mat¬ 
ter halos and subhalos have cuspy density distribu¬ 
tions that ar e we ll described b y NFW or Einasto pro¬ 
files ( see a lsolLovell et al .1120121) . iPontzen fe Governatol 
( 2012 ft andlShen et al.| (1 2014ft argue that baryonic physics 
( Governato et all 12012 b via a process similar to vio¬ 
lent relaxation, can lead to changes in the potential, 
specificall y, a more shallow potential and hence a core 
(see also iBrook et al.1 1201 2t iZolotov et al.112012ft . Al¬ 


ternatively, self-interacting dark matter can produce 


large cores ( 

Vogelsberger et al.1 20121 

Rocha et all 2013fc 

Zavala et al 

20131; iVoffclsberger et al. 

loil. 


4.2. Cusp-core transition 


Real systems like dwarf spheroidals may have shal¬ 
lower central logarithmic slopes and larger central re¬ 
gions of near constant density than pure dark matter 
halos. How does one make the central potential less 
deep than found in numerical simulations? T he change is 
likely brought about by baryonic physics dNavarro et al.1 
il996f : Governato et all 12012! : lOnorbe et al.l I2015T) which 
may induce a transformatio n from a typical —1 loga rith- 
mic slope to a cored halo dPeharrubia et al]l2012ft . In 
our picture this would correspond to a change in the 
central potential, through some energy injection pro¬ 
cess. For DARKexp to be valid, this energy injection 
would have to be collisionless to preserve the assump¬ 
tions behind the DARKexp derivation. It is interesting 
to note that the ba ryonic fe edback process proposed by 
IPontzen &; Governatol (|2012f) is not only collisionless but 
also dissipationless. 

We here illustrate the cusp-core transition as a trans¬ 
formation between two DARKexp halos. We can use 
DARKexp as a starting point for modeling baryonic ef¬ 
fects if they are mostly non-collisional, in which case 
DARKexp is also a good descriptor of the final state. 
Note, however, that we do not intend to model the de¬ 
tails of the transiti on in terms of detailed physics (e.g., 
lOnorbe et al.1 12015 1 such as the mass accretion history 
or star-formation history of the system, the supernova 
energy deposition efficiency, or the stellar initial mass 
function. 

Cusps or cores of dark matter halos are governed by 
the normalized dimensionless depth of the central poten¬ 
tial well, to the extent that their potentials are shaped 
by collective relaxation. How different are halos with 
different normalized central potentials, apart from their 
central logarithmic density slopes? We obtain analytical 
expressions for their total potential energies in Appendix 
A. We can also express the total potential energy of a 
halo as 


w = -c 


GM 2 

Vh 


( 8 ) 


with « 0.4 (lSpitzciifl969ft (see Figure (U). Here M is 
the total mass (which is finite for DARKexp and the 7 - 
models) and rh is the half-mass radius. 

For a 7 -model, the halo binding energy relative to an 
NFW cusp model, i.e., 7 CU s P = 1, is \AW/2W core \ = 



Y 



Ycore 

Figure 4. Left: The prefactor £ in Equation JHJ) for DARKexp 
(blue, solid) and the 7 -models (red, dashed). Right: Relative en¬ 
ergy required to transform an NFW cusp (logarithmic slope —1) 
into a core with central observed logarithmic slope of 7 CO re- We 
plot |AW/2Wcore| for DARKexp (blue , soli d) a nd the 7 -models 
(red, dashed), given by Equations dA4D and dA 6 D , respectively. In 
both panels, 7 for DARKexp was obtained from Equation ©. 

1/3(1 — 7 core) (Equation IA6I) . This relation is plot¬ 
ted in Figure [4] along with the corresponding relation 
based on DARKexp. The binding energy of a cored halo 
(7core = 0) is ~ 30% less than that of a cuspy halo 
(bcusp = 1 ), relative to the binding energy of the cored 
halo. For 7 core = 0.5, 7 cusp = 1.5 the difference is ~ 25%. 
These differences are significant, i.e., halos do not have 
universal structures or relative binding energies if they 
have different central cusp strengths. 

For a dwarf galaxy with virial mass M = 3 x 10 9 Mm, 
rh = 10 kpc, the total potential energy is W ft 3 x 10 5 ° 
erg. The cusp-core transition therefore requires ~ 10 55 
erg or ~ 10 4 supernovae (for a supernova energy input 
of ~ 10 51 erg) for 7 core = 0 (higher values of jeore re¬ 
quire less energy input). This corresponds to a total 
stellar mass produced of ~ 10 6 M 0 , with an uncer¬ 
tainty of about a factor of 2 depending on the initial 
m ass function. Thes e estim ates are consistent with those 
of iPenarrubia et~ahl (12012ft and iGarrison-Kimmel et al.l 
(12013ft . Because supernovae do not transfer all their ki¬ 
netic energy into the dark matter, the energy injection 
is less efficient and more supernovae are needed. On the 
other hand, the energy input needed is diminished when 
considering the actual star forma tion and mass assem¬ 
bly history of dwa rf spheroidals (jAmorisco et al.ll20l4 
iMadau et al.ll20 fdT) . Finally, observations do not strictly 
require 7 = 0 , further reducing the energy requirements. 
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As discussed in Section 3.2, dark matter halos of clus¬ 
ters of galaxies also appear to have c entral l ogari thmic 
slopes that are more shallow than —1 (ISand et al.l 12004 


IZitrin & Broadhurstlil2009t 

Newman et al. 

2013b 

laf). The 

effect seems to be generic 

Newman et al. 

2013a 

), giving 


central logarithmic slopes in the range 0.2-0. 8 . While 
some of this may be due to substructure, line of sight 
effects, or elongated structures (e.g., due to merging), 
it does appear that central profiles more shallow than 
NFW are observed in clusters of galaxies. However, for 
a massive cluster of galaxies with M = lx 10 15 M@ and 
rh = 500 kpc, the energy input required is ~ 10 6 ^ erg, 
corresponding to an unrealistic total stellar mass of 10 15 
Mg, if supernovae were to provide the energy input. It is 
therefore unlikely that supernova feedback is the origin 
of cores in cl usters of galaxies. It is possible that dynam¬ 
ical friction (lEl-Zant et al.ll200llf o f infalling satellites or 
black holes, active galactic nuclei dMartizzi et al.l [201^1 
or cluster me rgers dNewman et al.ll2013afl (although see 
IDehnenl I2005T) may provide the energy injection needed 
to erase the central cusp. Some clusters may also be born 
with a shallow potential, although it is noteworthy that a 
fair fraction of simulated cluster dark halos have central 
logarithmic slopes steeper than — 1 . 

5. CONCLUSIONS 

It has become standard practice to fit spherically aver¬ 
aged density profiles of dark matter halos, galaxies and 
clusters with phenomenological models, such as the NFW 
or its variants, or the Einasto a-model. We suggest that 
the reason these functions fit well is because they happen 
to resemble DARK exp, which is a theoretically derived 
model for collisionless self-gravitating systems, based on 
statistical mechanical arguments. 

Moreover, in this paper we have argued that non¬ 
universality of dark-matter halos is a natural expectation 
of DARKexp. If dark-matter halos would have logarith¬ 
mic slopes at small radii robustly resolved in simulations 
equal to — 1 , this would require some fine tuning of the 
model. The expectation of a range of properties appears 
to be borne out in simulations, and in observations of 
dwarf galaxies, and clusters of galaxies. In this picture, 
the range of central logarithmic slopes inferred can be 
seen as a reflection of a range in their central normalized 
potentials, or, equivalently, their total normalized bind¬ 
ing energies. This insight is useful when quantifying the 
extent to which halos differ as a result of their formation 
history and also can be used to analytically relate halos 
with different cusps to each other, such as required in the 
cusp-core transition. We do not suggest that the central 
potential is set up by some external process which then 
in turn causally sets the density profile. Rather, the cen¬ 
tral potential and the logarithmic slope of the density 
profiles are both signatures of an underlying formation 
process which determines the final equilibrium state. Fu¬ 
ture data on both dwarfs and clusters will be necessary 
to more fully examine the applicability of DARKexp to 
the whole range of systems. 

Because DARKexp does not have an analytical density 
profile, and if using tables of DARKexp p(r) is not an 
option, we have identified phenomenological fitting func¬ 
tions that match DARKexp reasonably well, and have 
presented relations between DARKexp’s <f>o and the best 


fitting 7 of the y-modcl, and a of the Einasto m odel. 
The resulting DARKexp dHiorth fe Wilhamsll2010H den¬ 
sity prof iles are well represented by Dehnen- Tremaine 7 - 
models (IDehnenl 119931 [Tremaine et al.lll994h . The main 
advantage of these models is that they allow for a range 
of values of the central cusp logarithmic slope and so can 
be used to model both the standard NFW-type halos 
with d log p{r)/d log r sa — 1 as well as more shallow (or 
steep) halos. 

As shown in th is paper (Figures CD and EH) and in - 
other recent wo r k (lOao et al]l2008 HSpringel et all[20081 
iVera-Ciro et al.l 120131: iLudlow et al.l 1 201 - 111 . cold dark 
matter numerical simulations by themselves generate a 
significant dispersion in a of simulated halos, reflecting 
different formation, accretion, or tidal stripping histories. 
DARKexp naturally accounts for such non-universality 
in profile shapes. In Figure [3] we demonstrate the non¬ 
universality of the pure dark matter halos. Both the 
shallower and the steeper systems are well fit by DARK¬ 
exp. While we have demonstrated a significant non¬ 
universality found in simulations of massive dark-matter 
halos, with a similar disp ersion noticabl e in simulations 
of lower mass halos (e.g.. ILudlow et aLll2013 ). it would 
be important to quantify the degree of non-universality 
and range of centr al potentials realized a s a function of 
halo mass (see also lDi Cintio et al.ll2014l ). 

In Section 4 we compare the potential energies asso¬ 
ciated with the core and cusp systems and discuss the 
required amount of energy and the possible injection 
mechanisms to bring about the transformation from a 
steeper pure dark matter system to a shallower observed 
system. The cusp-core transition cannot be simply mod¬ 
eled as a ‘one-size-fits-air transition from an NFW halo 
(or a DARKexp with </>o ~ 3 — 4) to a cored halo (or a 
DARKexp with <j>n < 2), because for both dwarf galaxies 
(jCollins et alJl2014 ) and clusters dNewman et al.ll2013af l 
there is a real dispersion in their observed properties. 
As highlighted in this paper the sa me is true for simu¬ 
lated massive da.rk-m atter halos (see ILudlow et al .1120131 
iDi Cintio et al.1120131 for other mass ranges). When 
modeling any transition of the central logarithmic slope 
properties as a result of some energy injection process, 
it is important to consider the significant dispersion in 
both the initial dark matter halos and those observed or 
affected by energy injection processes. 

Finally, we note that the ^>0 — 7 relation (Figure [l]) can 
be directly tested in numerical simulations, e.g., dissipa¬ 
tionless collapse or cosmological simulations. Halos with 
steeper profiles should have deeper central normalized 
potentials (possibly due to different formation times and 
accretion histories). For clusters of galaxies the <fi 0-7 
relation ca n be probed using gr avitational redshift mea¬ 
surements (jWoitak et aljl201lh of the central potential 
and gra vitational lensing measurements of the density 
profiles (lUmetsu et alJl201lt iBeraldo e Silva et alJl201& 
lUmetsu et al.ll 20 lblu 
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Cosmology Centre (DARK) is funded by the Danish Na¬ 
tional Research Foundation. L.L.R.W. would like to 
thank DARK for their hospitality. 
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APPENDIX 

A. POTENTIAL ENERGY OF CORED AND CUSPY DARKexp AND 7 -MODELS 
For DARKexp, N(e) = exp((jf > 0 — e) — 1, 


r<t> 0 

M(cf) 0 ) = / N(e)de = exp0 o - (</> 0 + 1) 
Jo 


and 


K + 2 W=- 


r<po 

/ eN{e)de 
Jo 


(Al) 

(A2) 

(A3) 


where K is the total kinetic energy and W is the total potential energy. Using the virial theorem, W = —2 K, 

^W = -M{fo) + ^-. 

Assuming unit-mass DARKexp models, the energy difference between a cusp and a core relative to the current core is 

( 2 - 


AW 


2 W r , 


M(4> o) 


cusp 


9 — 

For a 7 -model normalized to unit mass, p = ^r I ^r _7 (l + r) 7-4 , 


-1 • 


and 


W = —(10 — 47)' 


AW 1/5- 2 7c , 


2W CI 


2 \5 - 2j c 


-1 . 


(A4) 

(A5) 

(A 6 ) 


For the unit-mass 7 -models, the relation between central potential and inner logarithmic slope is 

07 = ( 2 - 7 ) _1 - 


The relevant radii are also directly related to 7 through 

r -2 = 1 - 7/2 

and the half-mass radius, 

r h = ( 2 4 /( 3 - 7 ) - I )" 1 

( Dehnen 1993 :1. Note that for -2 = 0.5 for 7 -models; for -2 ~ 0.5 also holds for a unit mass DARKexp, further 
demonstrating that the two models are similar. 


(A7) 

(AS) 

(A9) 


B. PHENOMENOLOGICAL FITS TO DARKexp 

Figure [5] shows a comparison between DARKexp and Einasto profiles. The correspondence is very good at small 
radii but poor at large radii. The poorer fit at larger radii is because the Einasto profile does not have a ‘built-in’ 
outer logarithmic slope of —4, except for values around the canonical a = 0.17. The profiles were matched in density, 
so the corresponding blue and green curves intersect each other at r a . The dashed curves in the same plot are the 
corresponding logarithmic density slopes; here too Einasto models track the respective DARKexp very well, although 
the blue and green dashed curves do not intersect each other at r a because the models can be matched either in density, 
or logarithmic density slope, not both. 

The left panel of Figure [ 6 ] shows DARKexp offo = 2,4, 6 as solid blue curves, and the matched 7 -models (Equation[3]) 
as red solid curves. These were matched in density. The 7 -models approximate DARKexp very well over nearly 4 
decades in radius. The right panel of Figure | 6 ] shows the same DARKexp density profiles with the 7 -models that were 
matched in logarithmic density slope, at r a ■ These too provide very good fits to DARKexp. 

Figure [7] shows the effect of varying /3 in the generalized Hernquist model (Equation [2j), between 0.5 and 2.0 (here 
we show model pairs matched only in density, not logarithmic slope). A value of f) = 0.7 appears to provide a slightly 
better fit than f3 = 1 but the improvement is not sufficiently significant to justify the introduction of another parameter. 
Note, however, that in the inner parts the correspondence with DARKexp is inferior to that of the Einasto profile, 
despite the fact that the generalized Hernquist model has an extra shape parameter. Overall, 7 -models appear to be 
useful representations of DARKexp. 
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log(r/r_ 2 ) 

Figure 5. Comparison of DARK exp and Einasto models. DARKexp models with </> o = 2,4, 6 (from bottom to top at small radii) are shown 
in blue (solid curves: density profile; dashed curves: logarithmic density slope). In green we show u-models matched at log(r a /r_ 2 ) = —1.5 
(dashed black lines) in density. 



log(r/r_ 2 ) log(r/r_ 2 ) 

Figure 6. Comparison of DARKexp and 7 -models (similar to Figure 2). In red we show 7 -models (Hernquist models with (3 = 1) matched 
at log(r a /r_ 2 ) = —1.5 lines in either density (left) or logarithmic slope (right). 


?- 
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cn 
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log(r/r 2 ) 


log(r/r 2 ) 



log(r/r 2 ) log(r/r 2 ) 


Figure 7. Effects of varying /3 between 0.5 and 2.0. Matching is done in density at log(r/r_ 2 ) = —1.5 like in Figure [ 6 ] (left). 
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